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Abstract 

The production of (anti-)deuteron and (anti-)^He nuclei in Pb-Pb collisions at = 2.76 TeV 

has been studied using the ALICE detector at the LHC. The spectra exhibit a significant hardening 
with increasing centrality. Combined blast-wave fits of several particles support the interpretation 
that this behavior is caused by an increase of radial flow. The integrated particle yields are discussed 
in the context of coalescence and thermal-statistical model expectations. The particle ratios, ^He/d 
and ^He/p, in Pb-Pb collisions are found to be in agreement with a common chemical freeze-out 
temperature of Tchem ~ 156 MeV. These ratios do not vary with centrality which is in agreement 
with the thermal-statistical model. In a coalescence approach, it excludes models in which nucleus 
production is proportional to the particle multiplicity and favors those in which it is proportional 
to the particle density instead. In addition, the observation of 31 anti-tritons in Pb-Pb collisions is 
reported. For comparison, the deuteron spectrum in pp collisions at = 1 TeV is also presented. 
While the p/TT ratio is similar in pp and Pb-Pb collisions, the d/p ratio in pp collisions is found to be 
lower by a factor of 2.2 than in Pb-Pb collisions. 
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1 Introduction 

Collisions of ultra-relativistic ions create suitable conditions for producing light (anti-)nuclei, because 
a high energy density is reached over a large volume. Under these conditions, hot and dense matter, 
which contains approximately equal numbers of quarks and anti-quarks at mid-rapidity, is produced 
for a short duration (a few 10^^^ s). The system cools down and undergoes a transition to a hadron 
gas. While the hadronic yields are fixed at the moment when the rate of inelastic collisions becomes 
negligible (chemical freeze-out), the transverse momentum distributions continue to change until also 
elastic interactions cease (kinetic freeze-out). 

The formation of (anti-)nuclei is very sensitive to the chemical freeze-out conditions as well as to the dy¬ 
namics of the emitting source. The production scenarios are typically discussed within two approaches: 
(i) The thermal-statistical approach has been very successful not only in describing the integrated yield 
of the hadrons but also of composite nuclei HKH . In this picture, the chemical freeze-out temperature 
Tchem (predicted around 160 MeV) acts as the key parameter. The strong sensitivity of the abundance of 
nuclei to the choice of Tchem is caused by their large mass m and the exponential dependence of the yield 
on the temperature given by exp(—ni/Tchem)- (ii) In the coalescence model, nuclei are formed by protons 
and neutrons which are nearby in phase space and exhibit similar velocities ||4l|5l. A quantitative de¬ 
scription of this process is typically based on the coalescence parameter Ba and has been applied to many 
collision systems at various energies |[6 UT3l . The binding energy of light nuclei is very small (around few 
MeV), so they can hardly remain intact during hadronic interactions, even if only quasi-elastic scatter¬ 
ing during the hadronic phase with temperatures between 100 MeV and 170 MeV is considered. When 
produced thermally at chemical freeze-out, they might break apart and be created again by hnal-state 
coalescence |[T4ll . It turns out that both, the thermal approach and the coalescence mechanism, give very 
similar predictions ifTSl . 

The production of light nuclei has attracted attention already at lower incident energies in heavy-ion 
collisions at the AGS, SPS, and RHIC lITbUTSl . A study of the dependence on .^^nn is of particular 
interest, because different production mechanisms might dominate at various energies, e.g. a formation 
via spectator fragmentation at lower energies or via coalescence/thermal mechanisms at higher ones. In 
all cases, an exponential drop in the yield was found with increasing mass of the nuclei |[T9l l20l . At 
RHIC and now at the LHC, matter with a high content of strange and of anti-quarks is created in heavy- 
ion collisions. This has led to the first observation of anti-alphas |[^ and of anti-hypertritons 1(2^ . Their 
yields at LHC energies were predicted based on thermal model estimates in l[T]|2l. 

In this paper, a detailed study of light (anti-)nuclei produced in the mid-rapidity region in Pb-Pb colli¬ 
sions at i/^nn = 2.76 TeV and a comparison with deuteron production in pp collisions at ^/s = 1 TeV 
using A Large Ion Collider Experiment (AEICE) ll23ll is presented. The paper is organized as follows: In 
Section]^ details of the analysis technique used to extract raw yields, acceptance and efficiency correc¬ 
tions of (anti-)deuterons and (anti-)^He are presented. The results are given in Sectionj^which starts with 
a comparison of the production of nuclei and anti-nuclei along with studies related to the hadronic inter¬ 
action of anti-nuclei with the detector material. Then, the transverse momentum spectra, px-integrated 
yields and average transverse momenta are shown. The observation of (anti-)tritons is also discussed 
in this section. In Section the results are discussed along with a description using a blast-wave ap¬ 
proach, and are compared with expectations from the thermal-statistical and coalescence models. The 
measurement of (anti-)alphas and (anti-)hypertritons will be shown in subsequent publications. 
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2 Experiment and data analysis 

2.1 The ALICE detector 

The results presented in this paper are obtained from the data collected by the ALICE experiment at the 
EHC. Its performance and the description of its various subsystems are discussed in detail in Refs, ll^ 
\2M . The AEICE detector has excellent particle identification capabilities. The main detectors used in 
this analysis are the Time Projection Chamber (TPC) ll25ll . the Time-Of-Elight detector (TOP) |[26l . and 
the Inner Tracking System (ITS) l[27l . All detectors are positioned in a solenoidal magnetic field of 
B = 0.5 T. As the main tracking device, the TPC provides full azimuthal acceptance for tracks in the 
pseudo-rapidity region |t 7 | < 0.8. In addition, it provides particle identification via the measurement of 
the specific energy loss AElAx. It allows the identification of (anti-)^He over the entire momentum range 
under study and the measurement is only limited by the available statistics. The velocity information 
from the TOP detector is in addition used to identify deuterons with transverse momenta (pj) above 
1.4 GeV/c and (anti-)tritons in the transverse momentum range of 0.6 GeV/c < px < 1-6 GeV/c. The 
detector provides a similar acceptance as the TPC and its total time resolution for tracks from Pb-Pb 
collisions corresponds to about 80 ps which is determined by the intrinsic time resolution of the detector 
and the accuracy of the start time measurement. By a combined analysis of TPC and TOP data, deuterons 
are identified up to 4.5 GeV/c in Pb-Pb collisions. In case of pp collisions, the less precisely determined 
start time leads to a time resolution of about 120 ps and the identihcation is limited to about 3 GeV/c. The 
precise space-point resolution in the six silicon layers of the ITS allows a precise separation of primary 
and secondary particles in the high track density region close to the primary vertex. 

2.2 Event and track selection 

Por this analysis, the data collected in the year 2010 are used. In total, the data sample consists of nearly 
14 million Pb-Pb collisions at i/^nn = 2.76 TeV and 380 million minimum-bias triggered events for pp 
collisions at ^/s = 7 TeV after off-line event selection. 

A pair of forward scintillator hodoscopes, the VO detectors (2.8 < rj <5.1 and -3.7 < rj < -1.7), 
measured the arrival time of particles with a resolution of 1 ns and was used for triggering purposes 
and for centrality determination of Pb-Pb collisions. In pp collisions, the data were collected using a 
minimum-bias trigger requiring at least one hit in either of the VO detectors or in the two innermost 
layers of the ITS (Silicon Pixel Detector, SPD). The trigger condition during the Pb-Pb data taking was 
changed with time to cope with the increasing luminosity delivered by the PHC. It was restricted offline 
to a homogenous condition, requiring at least two hits in the SPD and one hit in either of the VO detectors. 
This condition was shown to be fully efficient for the 90% most central events IITEII . A signal in the VO 
was required to he in a narrow time window (« 30 ns) around the nominal collision time in order to reject 
any contamination from beam-induced background. Only events with a reconstructed primary vertex 
position in the fiducial region |14| < 10 cm were considered in the analysis. The VO amplitude distribution 
was also used to determine the centrality of the heavy-ion collisions. It was htted with a Glauber Monte- 
Carlo model to compute the fraction of the hadronic cross section corresponding to a given range of VO 
amplitude. Based on those studies, the data were divided in several centrality percentiles, selecting on 
signal amplitudes measured in the VO ll28ll . The contamination from electromagnetic processes has been 
found to be negligible for the 80% most central events. 

In this analysis, the production of primary deuterons and ^He-nuclei as well as their respective anti¬ 
particles are measured at mid-rapidity. In order to provide optimal particle identification by reducing the 
difference between transverse and total momentum, the spectra are provided within a rapidity window 
of |y| < 0.5. In addition, only those tracks in the full tracking acceptance of |t]| <0.8 are selected. The 
extrapolation of the yield at low momenta, where the acceptance does not cover the full |y| < 0.5 region, 
is done by assuming a flat distribution in y and by determining drj/dy for each px-interval. Primary 
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Fig. 1: (Color online) Specific energy loss (dE/dx) vs. rigidity (momentum/charge) for TPC tracks from pp col¬ 
lisions at ^ = 1 TeV (top panel) and from 0-80% most central Pb-Pb collisions at .^/snn = 2.76 TeV (bottom 
panel). The solid lines represent a parametrization of the Bethe-Bloch curve. 
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Fig. 2: (Color online) Distribution of {m^ — m^) measured with the TOP detector for tracks with 2.6 GeV/c < pT < 
2.8 GeV/c from central Pb-Pb collisions showing the peak corresponding to the deuteron mass and the back¬ 
ground from mismatched tracks (black dotted line) which is subtracted to obtain the raw yields (see text for details). 


particles are defined as prompt particles produced in the collision including all decay products, except 
products from weak decays of light flavor hadrons and of muons. In order to select primary tracks of 
suitable quality, various track selection cuts are applied. At least 70 clusters in the TPC and two points in 
the ITS (out of which at least one in the SPD) are required. These selections guarantee a track momentum 
resolution of 2% in the relevant pT-range and a dElAx resolution of about 6%, as well as a determination 
of the Distance-of-Closest-Approach to the primary vertex in the plane perpendicular (DCA;(y) and par¬ 
allel (DCA^) to the beam axis with a resolution of better than 300 pm in the transverse direction ll24ll . 
Furthermore, it is required that the P^r TPC cluster is less than 4 and tracks of weak-decay products 
are rejected as they cannot originate from the tracks of primary nuclei. 

2.3 Particle identification 

Particle identification is mainly performed using the TPC |[25l . It is based on the measurement of the 
specific ionization energy deposit (dE/dx) of charged particles. Figure [T] shows the dE/dx versus rigidity 
(momentum/charge, p/z) of TPC tracks for pp collisions at y/s = 7 TeV (top panel) and for Pb-Pb colli¬ 
sions at -^Ann = 2.76 TeV (bottom panel). Nuclei and anti-nuclei like (anti-)deuterons, (anti-)tritons, and 
(anti-)^He are clearly identified over a wide range of momenta. The solid curves represent a parametriza- 
tion of the Bethe-Bloch function for the different particle species. In practice, it is required that the 
measured energy-loss signal of a track lies in a 3a window around the expected value for a given mass 
hypothesis. While this method provides a pure sample of ^He nuclei in the pj-vango between 2 GeV/c 
and 7 GeV/c, it is limited to about px < 1-4 GeV/c for deuterons. 

In order to extend the px-reach of the deuteron measurement, the TOF system is used above this mo¬ 
mentum in addition. Based on the measured flight time t, the mass m of a particle can be calculated 
as 


m 


2 




( 1 ) 


where the total momentum p and the track length L are determined with the tracking detectors. Fig¬ 
ure]^ shows the obtained Am^ distribution, where the deuteron mass square {m\) was subtracted, for a 
px-bin between 2.6 GeV/c and 2.8 GeV/c. For each px-bin, the Ani^ distribution is fitted with a Gaus- 
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Fig. 3: (Color online) Distribution of DCA;^, for deuterons (left) and anti-deuterons (right) in the transverse mo¬ 
mentum range 0.7 GeV/c < px < 1 -4 GeV/c for 0-80% most centtal Pb-Pb collisions at = 2.76 TeV demon¬ 
strating the influence of cuts in DCAj. on d and d. 


sian function with an exponential tail for the signal. Since the background mainly originates from two 
components, namely wrong associations of a track with a TOE cluster and the non-Gaussian tail of lower 
mass particles, it is described with a first order polynomial to which an exponential function is added. 
The same procedure for signal extraction and background subtraction is applied in the analysis of pp 
collisions. 

2.4 Background rejection 

Particles produced in the collisions might interact with the detector material and the beam pipe which 
leads to the production of secondary particles. The probability of anti-nucleus production from the 
interaction of primary particles with detector material is negligible, whereas the sample of nuclei may 
include primary as well as secondary particles originating from the material. This contamination is 
exponentially decreasing with increasing momentum. In addition, it is about five times larger in central 
compared to peripheral Pb-Pb or pp events because of the higher probability of a fake ITS hit assignment 
to secondary tracks. Most of the secondary particles from material have a large DCA to the primary 
vertex and hence this information is used to correct for the contamination. Eigure shows the DCAj^j, 
distribution for deuterons (left panel) and anti-deuterons (right panel) for Pb-Pb collisions at = 

2.76 TeV. The distributions are shown for two different IDCA^I cuts. As can be seen from the figure, 
a stricf jDCA^I cut of 1.0 cm cuts a large fraction of background for nuclei, but does not change the 
distribution for anti-nuclei. At sufficiently high momenta (above 1.4 GeV/c for deuterons and above 2 
GeV/c for ^He), the secondary and knock-out contamination caused by material is in this way reduced 
to a negligible level and the raw yield can be directly extracted. In order to extend the measurement of 
deuterons to lower momenta in Pb-Pb collisions, the DCAj^j, distribution for deuterons in each transverse 
momentum (pxl-interval was htted with the expected shapes (called “templates” in the following) as 
extracted from Monte-Carlo events. Eigure shows a typical example of this procedure for tracks with 
transverse momentum range 0.9 GeV/c < px < 10 GeV/c. One template for primary particles and one 
template for secondary particles from material are used. The characteristic shape of the template used for 
knock-out nuclei from material with its flat behavior at large DCAj^y allows a precise distinction between 
the two contributions. The significant peak at small |DCA;cj,| is caused by those knock-out nuclei to which 
a cluster in one of the SPD layers is wrongly associated. The obtained fraction of primary particles is 
then used to calculate the raw yield in the corresponding px-bin. The same technique is applied for 
background rejection and raw yield extraction of deuterons for pp collisions at ^/s = 7 TeV. 
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Fig. 4: (Color online) Distribution of DCA^y of identified deuterons in the transverse momentum range 0.9 GeV/c 
< Pt < LO GeV/c for central Pb-Pb collisions (.y/SNN = 2.76 TeV) along with the Monte-Carlo templates which 
are fitted to the data (see text for details). 


2.5 Efficiency and acceptance 


The final px-spectra of nuclei are obtained by correcting the raw spectra for tracking efficiency and accep¬ 
tance based on Monte-Carlo (MC) generated events. Standard event generators, such as PYTHIA ||29]| . 
PHOJET ||30ll . or HIJING ||3T]| do not include the production of (anti-)nuclei other than (anti-)protons 
and (anti-)neutrons. Therefore, nuclei are explicitly injected into underlying PYTHIA (in case of pp) 
and HIJING (in case of Pb-Pb) events with a flat momentum distribution. In the next step, the particles 
are propagated through the AEICE detector geometry with the GEANT3 transport code |[32l . GEANT3 
includes a basic description of the interaction of nuclei with the detector, however, this description is 
imperfect due to the limited data available on collisions of light nuclei with heavier materials. Due to the 
unknown interaction of anti-nuclei with material, these processes are not included for anti-nuclei heavier 
than anti-protons. In order to account for these effects, a full detector simulation with GEANT4 as a 
transport code 1331 |34l was used, hollowing the approach described in |[35l, the correction for interac¬ 
tion of (anti-)nuclei with the detector material from GEANT3 was scaled to match the expected values 
from GEANT4. An alternative implementation to correct for this effect and the relevant uncertainties 
related to these corrections are discussed in Section 3.1 The acceptance xefficiency is then obtained as 
the ratio of the number of particles detected by the detector to the number of generated particles within 
the relevant phase space. 


Eigure 1^ shows the acceptance xefficiency for deuterons (left) and ^He (middle) as a function of pj for 
Pb-Pb collisions at .^^nn = 2.76 TeV. In both cases, the rapid rise of the efficiency at low pj is determined 
by energy loss and multiple scattering processes of the incident particle with the detector material. The 
values reach a maximum when the energy loss becomes smaller and when the track curvature is still 
sufficiently large so that a track can cross the dead area between two TPC readout chambers in a relatively 
small distance such that the two track parts can still be connected. Eor straighter tracks at higher pj which 
cross the insensitive region between two chambers this distance is larger and the connection becomes 
more difficult. Thus a slight reduction of the efficiency is observed until a saturation value is reached. 
The figure also shows the lower efficiency values (open points) when in addition a deuteron track is 
matched to a hit in the TOP detector. The drop is mainly caused by the energy loss and multiple scattering 
in the material between the TPC and the TOP, by the TOP dead zones corresponding to other detectors 
or structures, and by the number of active TOP channels. The curves represent fits with the empirical 
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Fig. 5: Acceptance X efficiency as a function of transverse momentum (pj) for deuterons (left) and for ^He (middle) 
in Pb-Pb collisions at = 2.76 TeV, as well as for deuterons in pp collisions at y(s = 7 TeV (right panel). The 
curves represent a fit with the function presented in Eq. Q (see text for details). 



Fig. 6: (Color online) The average difference between the reconsttucted and the generated px is plotted as a 
function of the reconstructed px for simulated deuterons and ^He for Pb-Pb collisions at -y/iNN = 2.76 TeV. The 
lines represent a fit with the functional form as shown in Eq. ([^ (see text for details). 


functional form 


/(pt)=«o^^ ^+«3 Pt • (2) 

Here, aQ,a\,a 2 , and <33 are free parameters. Correcting the raw spectra with either the fit function or the 
actual histogram is found to result in negligible differences with respect to the total systematic error. 

Figure]^ (right) also shows acceptance xefficiency for the deuterons as a function of px for pp collisions 
at y/s = 7 TeV. The curve is a fit using the same functional form as used for the Pb-Pb collisions discussed 
above. The efficiency has a similar px-dependence as the one for Pb-Pb collisions at ^/^nn = 2.76 TeV. 
The observed differences are due to variations in the number of active detector components, mainly in 
the SPD, for the two data sets. 

2.6 Momentum correction 

Low-momentum particles lose a considerable amount of energy while traversing the detector material. 
The track reconstruction algorithm takes into account the Coulomb scattering and energy loss, assuming 









Production of light (anti-)nuclei in pp and Pb-Pb collisions at LHC energies 


ALICE Collaboration 


the pion mass for each particle. Therefore, a track-by-track correction for the energy loss of heavier 
particles (d/d and ^He/^He) is needed. This correction is obtained from MC simulations, in which the 
difference of the reconstructed and the generated transverse momentum is studied on a track-by-track 
basis. Figure shows the average px-difference as a function of the reconstructed track momentum 
for deuterons and ^He. The lines represent the empirical function 

/(pt)=A + b(^1 + ^^ , (3) 

where the free parameters A,B,C, and D are extracted from a fit. It can be seen that the correction 
becomes largest for the heaviest particles at low momenta. This reflects the typical l/jS^-behavior of 
the energy loss. The difference in transverse momentum is corrected on a track-by-track basis in the 
analysis. This energy loss correction has been applied both for pp and for Pb-Pb collisions. The same 
correction in rapidity has also been studied and found to result in negligible changes in the final spectra. 

2.7 Systematic uncertainties 

Individual contributions to the systematic error of the measurement are summarized in Table [T] and are 
discussed in detail in the following. The systematic uncertainty related to the identification of the nuclei 
is smaller in the px-region in which the energy loss in the TPC provides a clear separation compared to 
those in which the identification is mainly based on the TOE information. The error is of the order of 1% 
for deuterons at low momenta and for the full px-range studied for ^He-nuclei. In the TOE part (px > 
1.4 GeV/c) of the deuteron spectrum, the error is considerably larger due to the presence of background 
and has been estimated as 5% on the basis of different signal extraction methods: the raw yields obtained 
from the signal fit and from bin counting are compared. The estimates of the uncertainties related to the 
tracking and matching are based on a variation of the track cuts and are found to be less than 4% and 
independent of the particle species. In addition to this, a variation in the momentum correction leads to 
differences of similar magnitude at lower momenta and are added in quadrature. 

Contamination from secondaries originating from interactions of primary particles with the detector ma¬ 
terial dominates the systematic error at low transverse momenta, but it decreases exponentially towards 
higher momenta. These uncertainties are estimated by a variation of the fit range and templates. Their 
values amount to about 20% in the lowest px-bin for deuterons and for ^He in most central events. 
For all other centralities and transverse momentum regions, it is significantly lower. Feed down from 
weakly decaying hyper-nuclei is negligible for deuterons. The only relevant decay of the hyper-triton, 
d-|-p-|-7r^, results in a negligible contamination, because of the roughly 700 times smaller pro¬ 
duction cross section of the hyper-triton with respect to the deuteron mill. On the other hand, the decay 
—)■ ^He -I- 71^ contaminates the ^He-spectrum as these particles are produced with similar abundance. 
This background is conceptually similar to the feed down of A decays into the proton-spectrum 0^ 
though the relevant branching ratio in the case of (25%) llT7]l is assumed to be considerably lower 
than in the case of A (64%). A detailed MC study shows that only about 4-8% of all decaying into 
^He pass the track selection criteria of primary ^He. Therefore, the remaining contamination has not 
been subtracted and the uncertainty related to it was further investigated by a variation of the DCA;cj,-cut 
in data and a final error of about 5% is assigned. Uncertainties in the material budget have been studied 
by simulating events varying the amount of material by ±10%. This leads to variations in the efficiency 
of about 5% in the lowest px-bins. The hadronic interaction of nuclei with the detector material gives 
rise to an additional uncertainty of about 6% for deuteron and for ^He. The material between TPC and 
TOE needs to be considered only for the deuteron spectrum above px > 1.4 GeV/c and increases the un¬ 
certainty by additional 7%. The corresponding corrections for anti-nuclei are significantly larger and less 
precisely determined because of the missing knowledge of the relevant elastic and inelastic cross sec¬ 
tions. Details of the systematics originating from differences between the available models are discussed 
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Source 

d 

^He 


0.7 GeV/c 4 GeV/c 

2 GeV/c 8 GeV/c 

PID 

1 % 

5% 

1 % 

1 % 

Tracking and matching 

6 % 

4% 

6 % 

4% 

Secondaries material 

20 % 

1 % 

20 % 

1 % 

Secondaries weak decay 


negl. 

5% 


Material budget 

5% 

1 % 

3% 

1 % 

Hadronic interaction 


6 % 

6 % 



Table 1: Summary of the main contributions to the systematic uncertainties. See text for details. 


in the next section. 

In general, the individual contributions to the systematic error do not show a signihcant dependence on 
the event multiplicity. The only exception is given by the uncertainty of the correction for secondaries 
from material, which changes from about 20% in central to about 4% in peripheral Pb-Pb or pp collisions, 
respectively. All other contributions are found to be independent of event multiplicity. 


3 Results 

3.1 Anti-particle to particle ratios and hadronic interaction of anti-nuclei 

For a measurement of particle to anti-particle ratios, the correction of the hadronic interaction of the 
emitted particles with the detector material has to be precisely known. The relevant cross sections are 
only poorly measured for anti-nuclei heavier than p. The only available data for anti-deuterons from 
the U-70 Serpukhov accelerator ll^l^ are measured at relatively high momenta {p = 13.3 GeV/c and 
p = 25.0 GeV/c) and provide only a rough constraint. Two approaches are considered to model the cor¬ 
rection for hadronic interaction. Firstly, the anti-nuclei cross sections are approximated in a simplihed 
empirical model by a combination of the anti-proton (ap a) and anti-neutron (an^) cross sections. Fol¬ 
lowing the approach presented in BOll . the cross section for an anti-deuteron on a target material 
with mass number A is then e.g. given by 





(4) 


where the scaling factor K{A) is determined from the same procedure applied to the measured inelastic 
cross sections of nuclei and protons. Details of the method can be found in BOl . This approach is 
implemented as a modification to GEANT3. However, it does not account for elastic scattering processes 
and is therefore only used for the estimation of the systematic uncertainty. Secondly, the anti-nucleus- 
nucleus cross sections are determined in a more sophisticated model with Glauber calculations based on 
the well-measured total and elastic pp cross section Il34l . It is implemented in the GEANT4 software 
package 1331. 

The relevant correction factor for the anti-particle to particle ratio is given by the ratio of the efficiencies 
in which all effects cancel except of those related to the hadronic interaction with the detector material. 
The efficiency ratios for anti-deuterons and for ^He nuclei using the two models described above (mod¬ 
ified GEANT3 and GEANT4) are shown in Fig. [7] The applied correction factors are parameterized 
with the same function which was used for a similar study in l35l . The absorption correction is larger 
for tracks which are required to reach the TOF detector due to the additional material behind the TPC, 
mainly the support structure and the Transition Radiation Detector (TRD). In the following, results cor¬ 
rected with GEANT4 are presented. Based on the discrepancy between the two models, an uncertainty of 
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Fig. 7: (Color online) Ratio of anti-particle to particle efficiency based on GEANT4 and a modified version of 
GEANT3 including an empirical model to describe the hadronic interaction of anti-nuclei for (anti-)deuterons 
(left) and for (anti-)^He (right). The estimate of the systematic uncertainty for the hadronic interaction based on 
the difference between the two models is indicated by the blue band. 



Fig. 8: (Color online) Ratios of d and d as well as of ^He and ^He versus pi per nucleon for various centrality 
classes in Pb-Pb collisions at = 2.76 TeV. Boxes describe the systematic uncertainties, vertical lines the 
statistical ones. 


11 













































Production of light (anti-)nuclei in pp and Pb-Pb collisions at LHC energies 


ALICE Collaboration 


Anti-nuclei/nuclei 

Centrality 

Ratio 


0 -10% 

0.98 

±0.01 

±0.13 


10 -20% 

0.99 

±0.01 

± 

0.13 

d/d 

20-40% 

1.01 

±0.01 

± 

0.14 


40-60% 

1.02 

±0.01 

± 

0.16 


60-80% 

1.02 

±0.02 

± 

0.16 

^fk/^He 

0 -20% 

0.83 

±0.08 

± 

0.16 

20-80% 

1.03 

±0.14 

± 

0.18 


Table 2: Anti-particle to particle ratios for various centrality classes in Pb-Pb collisions at ^snn = 2.76 TeV. The 
first error represents the statistical error and the second one is the systematic error. See text for details. 


60% of the difference between the efficiency for particles and anti-particles is assumed for the absorption 
correction. It is indicated by the blue band in Fig. 

Applying this model-based correction to the data, leads to d/d and ^He/^He ratio shown in Fig. for 
various centrality bins in Pb-Pb collisions. Both ratios are consistent with unity and exhibit a constant 
behavior as a function of pj as well as of collision centrality. Since the same statements hold true for the 
p/p ratios HU . these observations are in agreement with expectations from the thermal-statistical and 
coalescence models El which predict a ratio of d/d = (p/p)^ and ^He/^He = (p/p)^. Table|^show the 
anti-particle to particle ratios for various centrality classes in Pb-Pb collisions at = 2.76 TeV. 

Ongoing studies on the hadronic interaction of anti-nuclei in the material between the TPC and TOF will 
allow to constrain the uncertainties of the currently purely model based corrections and to replace them 
with data driven ones. As the spectra for nuclei and anti-nuclei are consistent within the currently large 
uncertainties, only the spectra of nuclei are provided in the following. 

3.2 Spectra of nuclei 

The final spectra of deuterons obtained in Pb-Pb and pp collisions are shown in Fig. The statistical 
and systematic errors are shown separately as vertical lines and boxes, respectively. In pp collisions, the 
spectrum is normalized to the number of all inelastic collisions (Mnel) which includes a correction for 
trigger inefficiencies (see H2l [43l for details). It is fitted with the following function Il44l - i4^ that has 
been used for lighter particles 

1 <fN _ dN {n-\){n-2) f mx-mo V" 

InpjdpjAy dy In nC{nC+ mQ{n — l))\ nC ) 

with the fit parameters C, n, and the AN!Ay. The parameter mq corresponds to the mass of the particle 
under study (deuteron) at rest and mj = Pj to the transverse mass. As in the case of lighter 

particles, the function is found to describe the deuteron pj spectrum well in the measured range with 
a ;t;^/ndf of 0.26 . The fit function is used for the extrapolation to the unmeasured region at low and 
high transverse momenta (about 45% of the total yield) and a px-integrated yield of AN/Ay = (2.02 ± 
0.34(syst)) X 10^^ is obtained. 

While statistical errors are negligible, the systematic error is dominated by the uncertainty related to 
the extrapolation (13%) which is evaluated by a comparison of different fit functions HTI (Boltzmann, 
mj-exponential, px-exponential, Fermi-Dirac, Bose-Einstein). Based on the same extrapolation in the 
unmeasured region of the spectrum, a mean transverse momentum (px) of 1.10 ± 0.07 GeV/c is obtained. 

The final spectra of deuterons and ^He for Pb-Pb collisions at = 2.76 TeV are shown in Figs. 
and [T^ for various choices of the collision centrality. Again, the systematic and statistical errors are 
shown separately by boxes and vertical lines, respectively. The px distributions show a clear evolution. 
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(GeV/c) 


Fig. 9: (Color online) Efficiency and acceptance corrected deuteron spectra for Pb-Pb collisions at i/snn = 2.76 
TeV in various centrality classes and for inelastic pp collisions at y(s = 7 TeV. The dashed lines represent an 
individual fit with the BW function (Eq.|^ in the case of Pb-Pb spectra and with the function presented in Eq. Q 
in the case of the pp spectrum (see text for details). The boxes show systematic en'or and vertical lines show 
statistical error separately. 


becoming harder as the multiplicity increases. A similar behavior is observed for protons, which have 
been successfully described by models that incorporate a significant radial flow iHdl . 

The spectra obtained in Pb-Pb collisions are individually fitted with the blast-wave (BW) model for the 
determination of pj-integrated yields and (px)- This model BSl describes particle production properties 
by assuming that the particles are emitted thermally from an expanding source. The functional form of 
the model is given by 


I dN 

Pt dpx 


rdrmjlo 


px sinh p \ / mx cosh p \ 


Tkii 


'H 


Tkin 


) 


where the velocity profile p is described by 


( 6 ) 


p = tanh^' p = tanh^' (^j3s(r//?)”) . (7) 

Here Iq and Ki are the modified Bessel functions, r is the radial distance from the center of the fireball 
in the transverse plane, R is the radius of the fireball, j3 (r) is the transverse expansion velocity, j3s is 
the transverse expansion velocity at the surface, n is the exponent of the velocity profile, and Tkin is 
the kinetic freeze-out temperature. The free parameters in the fit are Tkin, j3s, n, and a normalization 
parameter. Here, we present two alternatives: fitting the two particles separately (Figs. [9] and [T0|) and 
simultaneously (Fig. [TT]). The extracted values of the kinetic freeze-out temperature and radial flow 
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Fig. 10: (Color online) ^He spectra for two centrality classes (0-20% and 20-80%) are shown for Pb-Pb collisions 
at -y/iiw = 2.76 TeV. The spectra are fitted individually with the BW function (dashed lines). The systematic and 
statistical errors are shown by boxes and vertical lines, respectively. 



Fig. 11: (Color online) The top panel shows the combined fit of deuteron and ^He spectra with the BW function 
for 0-20% centrality for Pb-Pb collisions at i/snn = 2.76 TeV. The systematic and statistical errors are shown by 
boxes and vertical line, respectively. The lower panel shows the deviation of the spectra from the BW fits. 
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Centrality 

(iS) 

Ain (MeV) 

n 

dN/dy 

{ pT ) (GeV/c) ;t:Vndf 


d(0-10%) 

0.630 ±0.003 

77 ±2 

0.75 ±0.05 

(9.82±0.04±1.58) x lO^^ 

2.12 ±0.00 ±0.09 

0.10 

d (10-20%) 

0.613 ±0.004 

96±2 

0.78 ±0.06 

(7.60±0.04±1.25) X lO^^ 

2.07±0.01±0.10 

0.07 

d (20-40%) 

0.572 ±0.004 

100 ±2 

0.96 ±0.07 

(4.76±0.02±0.82) x lO^^ 

1.92±0.00±0.11 

0.07 

d (40-60%) 

0.504 ±0.017 

124 ±7 

1.04 ±0.19 

(1.90±0.01±0.41) X IQ-^ 

1.63±0.01±0.09 

0.01 

d (60-80%) 

0.380±0.010 

108 ±3 

1.85 ±0.35 

(0.51 ±0.01 ±0.14) X 10-2 

1.29±0.01±0.14 

0.21 


^He (0-20%) 

^He (20-80%) 

0.572 ±0.006 

0.557 ±0.007 

101±61 

101 ±37 

1.02 ±0.02 

0.99 ±0.03 

(2.76 ±0.09 ±0.62) x lO^'* 
(5.09±0.24±1.36)x 10-^ 

2.83 ±0.05 ±0.45 

2.65 ±0.06 ±0.45 

0.49 

0.20 

d, -^He (0-20%) 

0.617 ±0.009 

83 ±22 

0.81 ±0.06 



0.32 


Table 3: Summary of extracted yields dN/dy and mean transverse momenta (pj) based on the BW individual fits 
performed on the spectra for Pb-Pb collisions at -y/iNN = 2.76 TeV. The hrst error on dN/dy and (px) represents the 
statistical error and the second one is the combination of systematic and extrapolation errors, added in quadrature. 
See text for details. 


velocity are discussed in more detail in the next section. The results of these fits are summarized in 
Table where the values of dN/dy and (px) are also reported. The dN/dy values are extracted by 
individually fitting the spectra with the BW model. The extrapolation to px = 0 introduces an additional 
error which is again evaluated by a comparison of different ht functions and amounts to about 6% for 
central and 13% for peripheral collisions for the deuteron yields. In the ^He case, it contributes about 
17% and 16% to the total systematic errors for the 0-20% and 20-80% centrality class, respectively. 



Fig. 12: (Color online) The production yield dN/dy of light nuclei as a function of the particle mass measured 
for 0-20% centrality class in Pb-Pb collisions at y/sNN = 2.76 TeV. The line represents a ht with an exponential 
function. 


Figure shows the production yields of p, d, and ^He measured in the centrality interval 0-20% in 
Pb-Pb collisions which follow an exponential decrease with the mass of the particle. The penalty factor, 
namely the reduction of the yield by adding one nucleon, is 307 ± 76. Such an exponential decrease has 
aheady been observed at lower incident energies starting from those provided by the AGS ifT^ [T^ |T9l 
M, yet with different slopes. 


The mean transverse momentum (px) values obtained for d and ^He are compared to those of light 
particle species for Pb-Pb collisions at ^/^nn = 2.76 TeV (from EH) in Fig. 13 The figure shows that 
the (px) increases with increasing mass of the particle. Such a behavior is expected if all the particles 
are emitted from a radially expanding source. 
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Fig. 13: (Color online) Mean ttansverse momentum {pj) as a function of particle mass for various centrality 
classes are shown for Pb-Pb collisions at i/snn = 2.76 TeV. 


3.3 Observation of (anti-)triton 



Fig. 14: (Color online) Scatter plot of — m?) measured with the TOP detector versus pj. Only those tracks 
are shown which pass the pre-selection done by applying a 3(7 cut on the TPC dE/dx. The px-region in which the 
candidates are identified on a track-by-track basis is shown as red box. 


The combined particle identification capability of the TPC and TOE also allows a track-by-track iden¬ 
tification of low momenta (0.6 GeV/c < pj < \.6 GeV/c) anti-tritons as illustrated in Eig. In this 


momentum region, the background from mismatched tracks is removed by the TPC particle identih- 
cation. The contamination is estimated based on a side-band study and found to be negligible below 
Pt < 1.6 GeV/c, but it increases rapidly for higher momenta so that signal and background cannot be 
distinguished anymore thus limiting the range available for the measurement. 


As can be seen, 31 anti-triton candidates are observed in the 0-80% centrality range. These numbers 
are consistent with expectations based on an extrapolation of the ^He-spectra to lower momenta taking 
into account the low reconstruction efficiency for anti-tritons in this momentum region (of about 11 % ± 
6 %). An observation of about 10 to 40 anti-tritons is expected based on this estimate, indicating similar 
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Fig. 15: (Color online) Blast-wave fit of K+, p, d, and ^He particles for 0-20% centrality for Pb-Pb collisions 
at -y/iNN = 2.76 TeV. Solid symbols denote the pj range of the spectra used for the hts, while the open symbols 
show the remaining part. The lower panels show the deviations of the measured spectra to the BW fits. 


production rates of anti-tritons and ^He nuclei. This comparison suffers from large uncertainties related 
to the absorption of anti-nuclei and energy loss in the detector material before the TPC at such low 
momenta. A similar measurement of tritons is unfeasible due to the large contamination from knock-out 
nuclei in this momentum region. 

4 Discussion 

4.1 Description of spectra via blast-wave fits 

Combined BW fits provide essential insight into the kinetic freeze-out conditions and allow quantitative 
comparisons between different collision systems and between measurements at different in terms 
of a hydrodynamic interpretation. In this section, a simultaneous fit to the 7i, K, p, d, and ^He spectra 
in the centrality range 0-20% using in addition data from ll4Tll4^ is discussed. Since the BW model is 
not expected to describe eventual hard contributions that may set in at higher pT, the fit ranges have been 
limited. For the light particles, they are taken as in BdlldOl (0.5-1 GeV/c, 0.2-1.5 GeV/c, 0.3-3 GeV/c 
for 71, K, and p, respectively). However, for d and ^He, the spectrum is fitted up to the pj value where 
the invariant yield reduces to 10% of the maximum available value of that spectrum. The exponent n 
of the velocity profile is left as a free parameter as in ll4T1l . In such an approach, all particle species are 
forced to decouple with the same parameters even though they feature different hadronic cross sections 
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with the medium. This is in particular relevant for multi-strange particles such as S and ifSOl . which 
are therefore not included in the fit. 


In Fig. 15 the results of a simultaneous fit to the five particle species are shown. The deviations of the 
spectra from the BW fit are shown in the lower parts of Fig. 15 The statistical errors are shown by 


vertical lines and the systematic errors are shown as shaded bands. Note that data points marked with 
open symbols are not included in the fit. The hardening of the spectra for central collisions is qualitatively 
well described by the combined BW fit with a collective radial flow velocity (j3) = 0.632 ± 0.01, a kinetic 
freeze-out temperature of Tkin = 113 ± 12 MeV, and n = 0.72 ± 0.03. The ;t;^/ndf value of the fit is 0.4. A 
comparison of these parameters to those obtained from a ht to n, K, and p BTI ((j3) = 0.644 ± 0.020, r^in 
= 97 ± 15 MeV, and n = 0.73 ± 0.11) reveals that the inclusion of nuclei leads to a slightly smaller value 
for (j8) and a slightly larger value for Tkin. This behavior is mainly driven by the strong anti-correlation 
of (j8) and in the blast-wave model: the slightly lower value of (j8) leads to a deviation of the fit 
from the proton spectrum which is then compensated by a higher Tjan. 


4.2 Comparison to thermal models 

Figure[^shows the d/p and the ^He/p ratios as a function of the average charged particle multiplicity per 
event. The proton yields are taken from 1441.14^ . The observed values of about 3.6 x 10^^ for the d/p ratio 
and about 9.3 x 10^® for the ^He/p ratio are in agreement with expectations from the thermal-statistical 
models mm. Similar values for d/p ratios are also observed by the PHENIX experiment for Au-Au 
collisions in 1211. Since at RHIC energies significant differences between nucleus and anti-nucleus 
production are present, for this plot the geometrical mean is used which in a thermal concept cancels 
the influence of the baryon chemical potential (/^b) *■ Within the achieved experimental precision, no 
dependence of these particle ratios on the event multiplicity is observed at RHIC and FHC energies. Also 
the p/p and the p/n ratios hardly vary with centrality 114111521 showing that Tchem and /Ib do not vary with 
centrality in high energy collisions. In a coalescence approach, the centrality independence disfavors 
implementations in which the nuclei production is proportional to the absolute proton multiplicity Il53l 
rather than the particle density. 


The comparison with thermal models is shown in more detail in Fig. 17 for the 0-10% centrality class. 
These calculations have been performed using the grand-canonical formulation of both THERMUS 1541 
and the GSI-Heidelberg model [HI. This approach is appropriate for the ratios shown here, as no strange 
quarks are involved. Details can be found in 111 |2l . These ratios are monotonically increasing with Tchem 
reflecting the dependence with exp(—Am/Tchem) where Am corresponds to the mass difference of the 
particles under study. 


The measured ratios of ^He/p and ^He/d are in agreement with a chemical freeze-out temperature in the 
range of 150 MeV to 165 MeV. No significant differences are observed between the THERMUS and 
GSI-Heidelberg model with respect to the production of light (anti-)nuclei. A fit to p, d, and ^He only 
gives Tchem = 156 ±4 MeV with a ;t;^/ndf of 0.4. This value can be compared to a fit including all 
measured light flavor hadrons which yields a temperature of about 156 MeV Il55l . 


At these temperatures, the weakly-bound deuteron and ^He can hardly survive. These nuclei might 
break up and might be regenerated. However, if this complex process of break-up and regeneration 
is governed by an overall isentropic expansion, the particle ratios are preserved ll56]l . Eventually, the 
yields of particles including weakly bound nuclei are therefore described in the thermal-statistical model. 
Other properties, e.g. spectral shapes and elliptic flow, exhibit the influence of the interactions during the 
hadronic phase. 

*In a thermal model, the yield ns of a baryon with energy £ in a medium of temperature T is proportional to exp(— ) 
while the yield of an anti-baryon n-g is proportional exp(—The geometric mean ^ngtig leads to a cancellation of the 

jUiJ. 
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Fig. 16: (Color online) d/p and ^He/p ratio in heavy-ion collisions as a function of event multiplicity. Within 
the uncertainties no significant variation with multiplicity is observed. The d/p and d/p results from the PHENIX 
Collaboration HtIIsTII are averaged as explained in the text. The lines represent fits with a constant to the ALICE 
data points. 



Fig. 17: (Color online) Particle ratios of nuclei as measured in 0-10% most central Pb-Pb collisions compared 
to the THERMUS ll54l model (solid lines) and the GSI-Heidelberg model IH (dashed lines) as a function of the 
chemical freeze-out temperature Tchem- The ^He yield is scaled to 0-10%. Horizontal error bars indicate the 
temperature range obtained by a projection of the total error of the ratio on the temperature axis. 
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The d/p ratio obtained in pp collisions is lower by a factor of 2.2 than in Pb-Pb collisions. Assuming 
thermal production not only in Pb-Pb, but also in pp collisions, this could indicate a lower freeze-out 
temperature in pp collisions. However, the p/Ti ratio does not show significant differences between pp 
and Pb-Pb collisions. Effects related to canonical suppression of strange particles can also be excluded 
because these ratios do not involve any strange quarks. Therefore, this observation must find anofher 
explanafion wifhin fhe framework of fhermal models or non-fhermal producfion mechanisms need fo be 
considered in small systems. Eurfher work in fhe fheorefical models is needed for a belter undersfanding 
of fhis effecf. 

4.3 Comparison with the coalescence model 

Eight nuclei have nucleons as constituents and are thus likely formed via coalescence of protons and 
neutrons which are near in space and have similar velocities. In this production mechanism, the spectral 
distribution of the composite nuclei is related to the one of the primordial nucleons via 


Ei 


d^Ni 

dp] 


= Ba 



( 8 ) 


assuming that protons and neutrons have the same momentum distribution. Ba is the coalescence param¬ 
eter for nuclei i with mass number A and a momentum of p, = A pp. 

Eigure[^ shows the obtained B 2 values for deuterons (left panel) and Bj, values for ^He (right panel) in 
several centrality bins for Pb-Pb collisions. The results are plotted versus the transverse momentum per 
nucleon. A clear decrease of B 2 and B 3 with increasing centrality is observed. In the coalescence picture, 
this behavior is explained by an increase in the source volume Vgff: the larger the distance between the 
protons and neutrons which are created in the collision, the less likely is that they coalesce. Alternatively, 
it can be understood on the basis of the approximately constant d/p and ^He/p-ratios as an increase of 
the overall proton multiplicity independent of the geometry of the collision. The argument can be best 
illustrated by assuming a constant value of B 2 and integrating Eq. ([^ over pj. The value of B 2 can then 
be calculated for a given ratio d/p and a given spectral shape /(px) (with /o°°/(pt) dpx = 1) of the proton 
spectrum as 


R =— ^ (Q\ 

where for a constant ratio of the deuteron dN^ /dy to proton dAp /dy yield, it is found that B 2 1 / (dAp/dy). 
As can be seen in Eig. the coalescence parameter also develops an increasing trend with transverse 
momentum for central collisions in contrast to expectations of the most simple coalescence models. The 
signihcance of this increase is further substantiated by the fact that the systematic errors between px- 
bins are to a large extent correlated. It can be qualitatively explained by position-momentum correlations 
which are caused by a radially expanding source lISTl . Eor quantitative comparisons, better theoretical 
calculations are needed which couple a coalescence model to a full space-time hydrodynamic descrip¬ 
tion of the fireball. Also in the discussion of the variation of the B 2 parameter as a function of collision 
energy, its strong dependence on centrality and px must be taken into account. It is observed that B 2 at a 
fixed momentum (px =1.3 GeV/c) for central collisions (0-20%) decreases rapidly from AGS energies to 
top SPS energy and then remains about the same up to RHIC lISTI . Our value of approximately 4 x 10^^ 
GeV^/c^ is only slightly lower than the measurement at RHIC (ps 6 x 10^"* GeV^/c^). 
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Fig. 18: (Color online) The coalescence parameters B 2 (left) and B 3 (right) as a function of the transverse momen¬ 
tum per nucleon for various centrality classes in Pb-Pb collisions at = 2.76 TeV. 


5 Conclusion 

In summary, the spectral distributions of deuterons in pp at y/s = 7 TeV and of deuterons and ^He in Pb- 
Pb collisions at = 2.76 TeV have been presented. In Pb-Pb collisions, the yields are decreasing 

by a factor of 307 ± 76 for each additional nucleon, the mean pT rises with mass and the combined 
blast-wave fit to n, K, p, d, and ^He gives a reasonable fit with (j 8 ) = 0.63 and Tkm around 115 MeV 
suggesting that the kinetic freeze-out conditions for nuclei are identical to those of the other light flavour 
hadrons. For anfi-frifons, a frack-by-frack idenfihcafion has been applied in fhe momenfum range 0.6 
GeV/c < pt < 1-6 GeV/c and fhe observafion of 31 anfi-frifons in Pb-Pb collisions af -^^nn = 2.76 TeV 
in fhe 0-80% cenfralify class is reporfed in Ibis paper. 

An imporfanf quesfion is whefher fhe nuclei produced in heavy-ion collisions are creafed af fhe chemical 
freeze-oul or af a lafer sfage via coalescence. One of fhe key observafions is fhe fad fhaf fhe d/p and 
^He/p rafios are consfanf as a funcfion of (dVch/d 77 ) 7 ^=o- Such a behavior is expecfed from a fhermal- 
sfafisfical inferprefafion, as if is found fhaf Tchem and Pb do nof vary wifh cenfralify in high energy 
collisions. Furthermore, a common freeze-oul femperalure of around 156 MeV for lighf (anfi-)nuclei 
and all ofher measured lighf flavor hadrons is obfained in a Ihermal-slafislical model. Fighf (anfi-)nuclei 
in Pb-Pb collisions Iherefore show fhe idenfical behavior as non-composife lighf flavor hadrons which 
are governed by a common chemical freeze-ouf and a subsequenf hydrodynamic expansion. 

The exfracfed coalescence paramefers B 2 and B 3 exhibif a significanf decrease wifh collision cenfralify 
and an increase wifh Iransverse momenfum which cannol be explained by coalescence models in fheir 
simples! form. On fhe ofher hand, faking info accounf fhe larger source volume for more cenfral collisions 
and fhe radial expansion of fhe emilfing sysfem, fhe producfion of lighf (anli-)nuclei in Pb-Pb collisions 
may sfill be compafible wifh fhe expecfalions from a coalescence space-fime descripfion. 

The measuremenfs of nuclei af FHC energies are shown fo follow frends observed from lower incidenf 
energies. Exfrapolafions and model predicfions based on fhe fhermal-sfafisfical or coalescence approach 
are, Iherefore, a solid ground for furfher sfudies, e.g. of hyper-nuclei and exofica. 
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